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Abstract 

In statistical mechanics, observables are usually related to local degrees of freedom such as the Q < 4 
distinct states of the Q-state Potts models or the heights of the restricted solid-on-solid models. In the con- 
tinuum scaling limit, these models are described by rational conformal field theories, namely the minimal 
models .M(p,p') for suitable p,p'. More generally, as in stochastic Loewner evolution (SLE K ), one can 
consider observables related to nonlocal degrees of freedom such as paths or boundaries of clusters. This 
leads to fractal dimensions or geometric exponents related to values of conformal dimensions not found 
among the finite sets of values allowed by the rational minimal models. Working in the context of a loop 
gas with loop fugacity |3 = — 2 cos we use Monte Carlo simulations to measure the fractal dimensions 
of various geometric objects such as paths and the generalizations of cluster mass, cluster hull, external 
perimeter and red bonds. Specializing to the case where the SLE parameter k = ^— is rational with 
p < p ', we argue that the geometric exponents are related to conformal dimensions found in the infinitely 
extended Kac tables of the logarithmic minimal models CM (p, p '). These theories describe lattice systems 
with nonlocal degrees of freedom. We present results for critical dense polymers CM[\ , 2), critical per- 
colation CM{2, 3), the logarithmic Ising model £M(3,4), the logarithmic tricritical Ising model CM (4,5) 
as well as CM{3, 5). Our results are compared with rigourous results from SLE K , with predictions from 
theoretical physics and with other numerical experiments. Throughout, we emphasize the relationships 
between SLE K , geometric exponents and the conformal dimensions of the underlying CFTs. 
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1 Introduction 

The numerical measurements of fractal dimensions and geometric exponents presented here aim at strength- 
ening the ties between three central paradigms describing two-dimensional phase transitions. The first, the 
lattice approach, which started with Onsager's calculation |1| of the spectrum of the transfer matrix of 
the Ising model and continued in the work of Baxter |2|, is the closest to a first-principle description of 
phase transitions. Rigorous results obtained from this approach are important, but remain far from the 
more physical continuum field-theoretic description used in solid-state or elementary particle physics. By 
limiting the description to the critical point, the second paradigm, conformal field theory (CFT) [3], has pro- 
vided such a full-fledged continuum field theory with exact calculations of certain correlation functions. In 
doing so, it enlarges the set of algebraic tools, already central in Onsager's solution, that describe critical 
behaviour. Even though CFT enables perturbation away from criticality, the approach to criticality is far 
from being understood. Moreover, its description of statistical models misses the probability framework 
that one would expect from models whose ab initio description uses the Boltzmann distribution. The third 
paradigm, Stochastic Loewner Evolution (SLE) l4ll5ll6ll7ll8ll9l lT0l , is indeed probabilistic in nature. In SLE, 
objects are of a measure-theoretic nature and are concrete enough to reproduce rigorously several predic- 
tions from the other two approaches, mostly those from CFT. Since rigorous derivations [11 1 give the fractal 
dimension of certain SLE paths confirming predictions from CFT, it seems natural to us that the study of 
fractal properties of critical geometric objects will improve the understanding of the ties between the three 
paradigms of lattice statistics, CFT and SLE. 
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Many geometric objects have been introduced to probe the critical behaviour of lattice models: the hull 
of a spin cluster, its mass, its external perimeter and its red bonds. In contrast to the usual correlation 
functions of local degrees of freedom, these geometric observables relate to nonlocal objects on the lattice. 
Theoretical predictions for various fractal dimensions and geometric scaling exponents were made in the 
late eighties llTZl . Many numerical measurements have been made, but only relatively recently were good 
measurements lTl3l [141 [TBI made of the fractal dimensions of these objects. These support the theoretical 
predictions but are limited to the Fortuin-Kasteleyn (FK) clusters [16| of Q-state Potts models \17\. 

The basic random variables in lattice spin models used for Monte Carlo (MC) simulations are the Potts 
spins sitting at the vertices of the lattice. Reasons to use Potts spin models for these simulations are nu- 
merous. First, spin variables are close to physical objects, readily interpretable as originating from an atom 
in a crystal. Second, the Swendsen-Wang algorithm |18| provides extremely quick upgrades in MC sim- 
ulations. Third, the thermodynamic limits of three of the Q-Potts models are archetypal lattice statistical 
models, namely, critical percolation (Q — > 1) with central charge c = 0, the Ising model (Q = 2) with c = A 
and the 3-state Potts model with c = 4 . 

The continuum scaling limits of the Q — 2,3 Potts models are described by rational CFTs. More explic- 
itly, they are described by the members M{3,4) and M{5, 6) of the series of minimal models M[y,y') HU 
20 1 with central charges 

c = c(R) = 1 - 6(1 - K)2 = 13-6(k+ -), k = ^>1, p,p' > 1 coprime (1.1) 

K x K' p 

related to Q by the mapping 

v/Q=-2cos?, Q<4. (1.2) 

K 

But now we are immediately presented with two problems. First, critical percolation, which formally is 
an allowed minimal model corresponding to p = 2, p ' = 3, is empty in the sense that it only contains the 
identity with conformal dimension A = 0. Second and more generally, the conformal dimensions related to 
the nonlocal geometric observables are not found |12| among the admissible conformal dimensions in the 
Kac tables for these rational CFTs. 

To circumvent these problems in considering fractal geometric objects, we follow II2T1I221 and move the 
lattice context away from the rational Q-state Potts or Restricted Solid-On-Solid (RSOS) models 120] and 
instead consider a loop gas with loop fugacity 

$=-2cos n . (1.3) 

K 

We believe there is good ground to measure the fractal dimensions of suitable geometric objects directly in 
the loop gas. In the loop gas, the basic random variables are the nonlocal degrees of freedom associated 
with segments of the loops. The continuum scaling limit of such a loop gas is not described by a rational 
CFT — it is described by the Coulomb gas f23l Moreover, if k is rational, which is the case of prime 
interest here, the loop gas is equivalent to the logarithmic minimal model CAi (p , p ' ) 1 24J and exhibits some 
remarkable properties. In the continuum scaling limit, it is thus described by a logarithmic CFT |25| with 
central charge given by the same expression as in the rational case ([TTTJ 

c = 1 - 6(1 - = 13-6(k+ -), k= 13 ->1, p,p' > 1 coprime. (1.4) 

K v K' p 

It is noted, though, that the logarithmic minimal model £.M(p,p') is defined also for p = 1. Critical dense 
polymers, with c = —2, now appears as the member £.M(1,2) and critical percolation as the member 
£jM(2, 3). For recent studies of critical percolation as a logarithmic CFT, we refer to [26J and references 
therein. As we will argue in many cases, the conformal dimensions related to the observed geometric 
exponents are found among the conformal dimensions in the infinitely extended Kac tables | |2.3) for these 
logarithmic CFTs. 

The basic probabilistic object of chordal SLE K is a path or curve with (at least) its starting point on the 



boundary of the domain within which it evolves. The evolution is described by 1 2.1 1 where we identify 



k =4k. 
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(1.5) 



Probabilists working with SLE use k = 6 for critical percolation while physicists working with CFT use 
i< = j. One can argue that studying a lattice loop gas is to probe directly a discretized version of SLE paths. 
Alternatively it has been proved for some values of k that the curves formed by loop segments of the loop 
gas coincide, in the continuum scaling limit, with SLE paths. Throughout this paper, we will emphasize the 
relationships between SLE, geometric exponents and the conformal dimensions of the underlying CFTs. As 
already indicated, for rational values of i< = p'/P/ we argue that the relevant underlying CFT is that of the 
logarithmic minimal model £A4(p,p'). 

The paper is organized as follows. The next section recalls the definition of the SLE process and sum- 
marizes the theoretical predictions for the fractal dimensions and geometric exponents of interest. Section|3] 
is devoted to the measurement of the fractal dimension of a path associated with a defect. This is the most 
natural discretization of an SLE path. Section [4] introduces natural extensions of classical geometric objects 
(cluster hull, cluster mass, external perimeter and red bonds) to the context of a loop gas. The measurement 
of these quantities by MC methods is also presented in this section. The lattice loop gas and the logarith- 
mic minimal models LM. (p , p ' ) are recalled in Appendix A. The various MC upgrade algorithms that we 
use are discussed in Appendix B. Finally, in Appendix C, we describe the initial thermalization of the MC 
process as well as checks on the statistical independence of our measurements. 

2 SLE and CFT 

2.1 SLE 

Two landmark results led to the development of SLE as a tool to probe critical phenomena. One of them is 
Smirnov's proof [27] that one observable of critical percolation on the triangular lattice has a scaling limit 
that is conformally invariant. The second is the introduction by Schramm [4] of a one-parameter family 
of conformally invariant measures on curves in the half -plane, denoted by SLE K . Schramm also showed 
that this family exhausts all such conformally invariant measures. Hence, the measure on interfaces of 
percolating clusters must fall in this family. Schramm showed that SLEg was the one. 

A simple way to understand this measure it to see the interfaces (or curves) as being grown out of a 
point on the real line, the boundary of the half -plane. During the growth, the interface may touch itself or 
the real line. Points belonging to the interface and points from which a curve cannot be drawn to infinity 
without crossing the interface form the SLE hull. 

Let us denote by Ct the complement of the hull at time t. Then Ct n H is a simply-connected open set 
and, by the Riemann mapping theorem, there exists an analytic function g t : C t D H — > H that is one-to-one. 
It is this map g t that satisfies the stochastic Loewner equation: 

2 

3 t g t z) = — — ?= — , t>0, withg (z)=z, (2.1) 

g t (z) - v K t>t 

where B t is a one-dimensional Brownian motion starting at 0. To determine uniquely g t , its behaviour 
around infinity is fixed by lim^oo (gt(z) — z) = 0. The interface y = Yt is then defined as the (continuous) 
path y : [0, oo) — > H with gt (yt) = >/icBt- Equation | |2.1| determines (somewhat implicitly) the probabilistic 
properties of yt in terms of those of the Brownian motion. 

One of the important corollaries of stochastic Loewner evolution is that the curve y generated by the 
process with parameter k has Hausdorff dimension equal to 

d^ = min(2,1 + £] (2.2) 

almost surely. This result, due to Beffara 1111 , together with the fact that the scaling limit of percolation 
is SLEg, proves one of the important predictions of theoretical physics, at least for one model on a very 
particular lattice. If one believes in universality, then the result should also apply to other lattices. Note 
that, if one sets k = in ( |2.1) , the measure degenerates and the interface yt tracing the imaginary axis has 
probability one. The fractal dimension of y is then 1 . For k larger or equal to 8, the curve is space filling and 
has dimension 2. 
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2.2 Logarithmic minimal models and CFT 

The logarithmic minimal models CM(y,y r ) [24] are reviewed briefly in Appendix |A| The degrees of free- 
dom of these theories are nonlocal objects associated with segments of the loops. In the continuum scaling 



limit, these models are described by logarithmic CFTs with central charges 1 1.4 1. Conformal dimensions of 
these CFTs can be organized into infinitely extended Kac tables 

A r ,^ (Rr - S)2 4 T (R - 1)2 , r.seN. (2.3) 
4k 

There is no claim here, however, that these values exhaust all possible values of the conformal dimensions 
for these theories. Indeed, there is some evidence l[T2l that conformal dimensions corresponding to half- 
integer values of r and s are allowed. 

There is a duality in these theories, mimicked in SLE K and described by 

16 1 

pnp , or kh — , or kh 7 . (2.4) 

K K 

Under this duality, the central charges and conformal dimensions verify 

c(R)=c(l), A TiS (I) =A s , r (R). (2.5) 



Here we only consider the case p < p ' with R = ^- > 1 . From the lattice, the dual models with p > p ' and 

R < 1 are expected to come from the dilute logarithmic minimal models with central charge c(k) = c (^) 
but we do not consider these models here. 



2.3 Relation between geometric exponents and conformal dimensions 

An important challenge is to relate the values of various geometric exponents or fractal dimensions d of 
observables to the conformal dimensions A of correlators 

(cp(O)cpM) - |x|- 2A (2.6) 

in the underlying CFTs, for example the critical Q-state Potts models. For Q integer, these are rational CFTs 
with relatively few conformal dimensions and sometimes one is led to conformal dimensions not found 
among the admissible values for these rational theories. Some observables in the logarithmic minimal 
models CM. (p , p ' ) should correspond to observables in the Potts models related to FK clusters. In addition, 
however, the logarithmic minimal models £.M(p,p') allow us to explain many observables not explained 
by the rational models. Consider a geometric observable O, measured on FK clusters of linear scale or 
radius R and described by correlators of the field cp>. An example is C, the cluster mass (total number of 
sites) of a percolating cluster, or H, its hull. A heuristic scaling argument in two dimensions then gives 



O(R) 



d 2 x(cM0)cMx)) 



d 2 x|*r 2A ~R 2 ~ 2A ~R do (2.7) 



where do stands for the fractal dimension of the observable O. In this way, such observables are associated 
with particular conformal weights in the Kac table. The renormalization group or geometric exponent y is 
related to A by 

y=2-2A = d . (2.8) 

In Section |3J we verify, through MC sim ulations, that precisely the fractal dimensions given in | |2.2| 
are found in the continuum scaling limit of the logarithmic minimal models £.M(p,p') with boundary 
conditions imposing a single (chordal) defect propagating through the system. We find it natural to identify 
this defect with (a discretized version of) the SLE path. We also note that for all p , p ', the dimension d^Jj, 
given above corresponds to 

d^ t E h=2-2A p , p ,_ 1 = W _, (2.9) 
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where A PiP /_i is a conformal dimension 1 2.3 1 of the logarithmic minimal model that falls outside the ratio- 
nal Kac table. 



In Section 4.1 we recall the definitions of the hull, cluster mass, external perimeter and red bonds for 
the context of the FK clusters of the Q -state Potts models and show how these definitions can be extended 
to logarithmic minimal models £.M(p,p'). The fractal dimensions (cluster mass) of FK clusters are given 
by (see IfTOl fT5 1 and references therein) 

— a a — 2 16 4k ' a C — 1 + 8 + 2k 

with the identification 

C = {Cluster mass} <-> A a = A(-p±i)/2,p'/2- (2.11) 

Notice that, depending on the parities of p and p', the Kac labels here can take half-integer values as 
foreshadowed after [ |2.3} . We also have the following associations (see ||T0ll28] and references therein) 

H = {Hull} <-> A P)P /±i , EP = {External Perimeter} <-> A p±1 , pS RB = {Red Bonds} <-} A VtV >±j (2.12) 

with conformal and fractal dimensions 



A FK 


= ^p,p'±1 


i 

~~ 2 


K 

4 » 


Q H 




(2.13) 


A FK 
Zi Ep 


= A p ±i >P ' 


1 

— 2 


1 

4k ' 




= 1+A 


(2.14) 


A FK 
^RB 


= A p , p '±2 


1 

~~ 2 


K i 3 

4 ' 4k' 


Q RB 




(2.15) 



We summarize in Table [T] the values of the fractal dimensions of these observables for the logarithmic 
minimal models CM{p,p') that will be studied. In Section |4j we verify, through MC simulations, that 
precisely these fractal dimensions are found by looking at suitably defined nonlocal observables of these 
logarithmic minimal models CM (p , p ' ) . We note that the fractal dimension of the cluster hull dn coincides 
with the fractal dimension of the SLE path dj^. 



3 Fractal dimension of a defect 

3.1 Definition and measurements 

For the simulations, we use the logarithmic minimal model CM{y,y') r a loop gas introduced in Il24l . Some 
very rudimentary properties of these models are recalled in Appendix [A] When the parameter R charac- 
terizing the model is chosen to describe one of the Q -Potts models, the distribution of the loops coincides 
with that of the FK contours. We refer to an open curve in a loop configuration, that is, one starting and 
ending on the boundary of the domain considered, as a defect. Because of the parallel between loops in the 
logarithmic minimal models and interfaces of FK graphs, it is natural to consider a defect as a (discretized 
version of a) chordal SLE process. The goal of the present section is to measure the fractal dimension of such 
a defect for models corresponding to various values of R of the SLE process. In the case of critical dense 
polymers CM (1,2) where R = 2, the loop fugacity vanishes, (3=0, hence excluding all configurations with 
at least one closed loop. In this model, the defect fills the space and the agreement with d^f h — 2 is trivial. 
This model is not studied further in this section. 

There are various definitions of dimension for objects of fractal nature. The Hausdorff dimension is 
widely used in analysis, but it is customary to use the Minkowski definition for simulations. Let D be a set 
Dcl d and let N ( e ) be the number of boxes of side length e required to cover D in an evenly-spaced grid. 
The Minkowski dimension of D is then given by 

do = lim (3.1) 
e->olog(1/ej 

if the limit exists. There are sets whose Hausdorff and Minkowski dimensions differ, but this subtlety is 
somewhat irrelevant here since we are to approximate the measurement of fractal objects by studying^infte 
geometric objects. 
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model 


K 


K 


C 




d H 


d c 




dRB 


CM( V ,V') 


4p' 
P 


vL 

V 


13-6R-| 


-2cosf 


1+f 




2k 


1 + 1 - ^ 

2 2k 


CM[1,2) 
CM(3,5) 
£M[2,3) 
£M{3,4) 
£M(4,5) 


8 

20 
3 

6 

16 

3 

5 


2 

5 
3 

3 
2 

4 
3 

5 
4 


-2 

3 
5 



1 

2 

7 

10 



1 

i(V5 + 1) 


2 
1 1 

6 

7 
4 

5 
3 

13 

8 


2 

77 
40 

91 
48 

15 

8 

299 
160 


5 
4 

13 
10 

4 
3 

1 1 

8 

7 
5 


5 
4 

14 
15 

3 
4 

13 

24 

17 

40 



Table 1: The models, their parameters k, R, c, |3 and fractal dimensions. The models include critical dense 
polymers CM. (1,2), critical percolation CM. (2,3), the logarithmic Ising model CM. (3,4) and the logarith- 
mic tricritical Ising model CM(4,5). The notations H,C,EP, RB stand for Hull, Cluster mass, External 
Perimeter and Red Bo nds . The fractal dimension of the hull coincides with the fractal dimension of the SLE 
path as given in |2.2|. The fractal dimensions of the cluster, external perimeter and red bonds coincide 
with those of the FK clusters in the associated Q -state Potts models via ( |1.2) . 




Figure 1: A configuration onaHxV = 9x8 lattice. The defect visits Na = 104 quarter-circles, including 
the ones on the boundary. 

The objective here is to examine, using simulations, whether the fractal dimension of a defect do indeed 
corresponds to the fractal dimension of the SLE path cffi^ 

d D I d*&. (3.2) 

However, since we work on a finite lattice, say of size H x V, we want to consider a discretization d^ x v of 
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the fractal dimension defined on the finite lattice yet approximating do in the infinite limit. One may think 
of the object D (the defect) as being drawn on a rectangle whose size and aspect ratio are fixed; the increase 
in H and V, while keeping their ratio fixed, then amounts to the limiting process. Therefore, if the size of 
the rectangle is fixed, H x V is related to the side length e of a box asHx V-^or 

log ( 1 / e ) ~ - log (# of boxes ) . (3.3) 

In fact, as each box can be crossed twice, one should instead count the number of quarter-circles, that is, the 
number of half -boxes. We therefore consider 

HxV logN A (D) log Na(D) f , 4 . 

d ° log(2HxVP log(R) (d) 

where Na(D) is the number of quarter-circles the defect D visits and where we have introduced some 
linear scale R. Up to a linear factor, R should be approximately (2H X V) i . If the defect crosses the same 
box twice, this counts for 2 in N&(D). Figure [l] depicts the setting that we use for measurements of cIq XV . 
The lattice is chosen to be almost square as the number H of boxes in a row is equal to V + 1 where V is the 
number of boxes in a column. We take V to be even. This allows for all loops on vertical boundaries to be 
closed by half -circles, as depicted in the figure. Because H = V + 1 , all loops but one reaching a horizontal 
boundary can also be closed by half -circles. If we choose the pair of loose ends at the top right and bottom 
left, respectively, the defect will cross the lattice and be a macroscopic object. It is for this geometric object 
that Na is computed for each configuration of the sample. To probe the asymptotic limit do, we will have 
to repeat these measurements on several lattices, starting with very small ones H x V = 3 x 2 and proceeding 
to 51 3 x 51 2 or 1 025 x 1 024, depending on the model. 

Here we present the results. 



The technical details of the statistical analysis are given in Section 3.2 



We study four models: critical percolation CM. {2, 3) and the logarithmic 



sing model CM (3, 4) but also 



two which do not have Potts model cousins, namely the model CM (3, 5) with c = — f and k = and 
the logarithmic tricritical Ising model CM(4,5) with c = and k = 5. The results for the latter models 
demonstrate that fractal dimensions can be measured also for models not having Potts model cousins. For 
all four models, measurements are made of d^ x v for several lattice sizes in two different ways, as described 
in Section 3.2 The results are plotted in Figure|2] Fits are made for both sets, discarding the smallest lattices. 



The full curves go from the smallest to the largest lattices used for the fit. The fitted curves are extended 
with dashed lines to the rejected measurements. Numerical results appear in Table [2] The datum 1 .835|1, 
for example, means that the 95%-confidence interval is of 1 unit of the last digit of 1 .835. We stress that this 
is the interval due to statistical uncertainty of the measurements. The asymptotic behaviour is very good 
for the models CM[3,5), CM(2,3) and C M (3 , 4) . There is a slight departure from the predicted theoretical 
value in the case of the logarithmic tricritical Ising model CM(4,5). One notices that the 95%-confidence 
interval does not always contain the theoretical prediction. The analysis in |3.2.2 however, will show how 
delicate the implied extrapolation is. 

The setting for our simulations is quite different from those used in previous measurements IIT31 ITU . 
Here are the main characteristics of our setting. 

• The linear size R of the objects is determined by the lattice size HxV. 

• All N^'s considered for the averages characterize objects spanning the whole lattice. These objects 
might not be the largest among all loops but, with high probability, will be among the largest. In a 
sense, they describe only percolating objects. 

• Also with high probability, the objects will bounce more than once on the boundary. They will there- 
fore probe both the bulk and the boundary and hence carry information about the conditions of both 
regions. 

It is interesting to recall the main features of the other simulations ||T31 ITSl alluded to above and contrast 
them with those of ours. 
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1024 256 64 
512 128 



32 



16 



"~3 



X LM(3,5) 



O LM(2,3) 



□ LM( 3,4) 
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Figure 2: The data {d^*^} (lower curve) and {dQ^ v } (upper curve) for the four models studied. The curves 
are the fits obtained as described in Section 3.2.2 to which end only lattices with H x V > 9 x 8 are used. 
Here and throughout the paper, the error bars are much smaller than the symbols used to depict the data. 



model 


K 


jtheo _ I i k 
U D — 1 + g 


dD,i 


dD,2 


CM(3,5) 


20 
3 


^ - 1.833... 


1.835|1 


1.835|1 


£M(2,3] 


6 


| = 1.75 


1 .75059162 


1.75133|58 


CM{3,4) 


16 

3 


1.667... 


1.665|1 


1.666|1 


£M(4,5] 


5 


^ = 1 .625 


1.618|1 


1.620|1 



Table 2: Values do,i and do, 2 with 95%-confidence intervals. 



Asikainen et al [13] perform measurements on a single square lattice of linear size 4096. In fact, they use 
two sets of measurements as they probe both open and periodic boundary conditions. They consider all FK 
graphs contained in a configuration. For each of these graphs, they calculate both its mass (the equivalent 
of Na above) and its gyration radius. The masses of all FK graphs of all configurations having their radius 
in a given window [R|, Rt+i] will contribute to this linear size bin. The characteristics of their setting are 
the following. 
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• The linear size R depends on the FK graph being studied, not on the lattice. 

• All objects of every size in a configuration are considered. 

• Smaller FK graphs may completely avoid the boundary and probe only the bulk. In the open bound- 
ary case, larger ones are likely to explore regions close and far from the boundary. 

While our approach is close in spirit to the growth process described by SLE, one can say that the study by 
Asikainen et al concerns global statistical features of FK graphs. 

The simulations carried out by Janke and Schakel 1 14J are limited to the Ising model but aimed at clean 
measurements. Like us, they perform measurements on several lattices. These are all square and their 
linear sizes range from 8 to 512. They measure both spin and FK clusters and, for each type of cluster, 
they count the number of sites occupied by the objects. This is a natural thing to do for spin clusters, but 
perhaps a little less natural for FK clusters. Contrarily to the Asikainen collaboration and us, they do not 
tackle directly the fractal dimensions of the objects they study. Instead, they extract the fractal dimensions 
from the percolation strength (the fraction of sites in the largest cluster) and the average cluster size x- 
This trick allows them to obtain very clean values for the fractal dimensions they measure from relatively 
small lattices. In summary, their setting is characterized by the following. 

• The properties (Poo and x) are studied as functions of the linear size of the lattices. 

• Pqq probes one particular object, similar, yet not equal, to the defect we are concentrating on. The 
observable x is an average over all clusters. 

• They use only periodic boundary conditions and therefore all their objects are in the bulk. 

As one can see, the three experiments (Asikainen et al's, Janke and Schakel's, and ours) are quite differ- 
ent. The fact that they nevertheless give very similar results is an indication that the physical observables 
considered are rather robust, with variations in boundary conditions or even definitions not causing major 
differences between measured fractal dimensions. There is no doubt, however, that the speed of conver- 
gence to asymptotic behaviour obtained by Janke and Schakel is impressive. We note that an earlier paper, 
that of Fortunato [29 1, also reports the cluster mass for the Ising and the 3-Potts models for the spin clusters. 
For the Ising model, Fortunato reproduces the predicted value ( ^£ that appears in Table [TJ to four signif- 
icant digits. For his measurements, he uses only the largest cluster in a configuration; this is not exactly 
what we do (ours always "percolate" from top to bottom but they might not be the largest) but his setting 
is nevertheless the one closest to ours. 



3.2 Technical issues 

3.2.1 Probability distribution of d^ xV and error on measurements 

The most obvious way of estimating^ d^ x v is 

-T£v 1 v logN A (DHxv ] 

d D,l =^2_ l ogR HxV (3-5) 

where the sum runs over the sample S with sample size |S|, R HxV is the total number of quarter-circles 
accessible to the defect (including those on the boundary), while N a ( D ^ x v ) is the actual number of quarter- 
circles visited by the defect in the configuration D[ J lxV . For a defect entering at one corner of a rectangle 
and exiting near the opposite one, the number N a (Dj^ x v ) is never zero and this definition works well. For 
other geometric objects to be studied in Section |4j however, a naive extension of this definition fails. For 
example, there are configurations D J, 1 x v whose number of red bonds is zero and the logarithm cannot be 
taken for this D[ J lxV . An alternative is 

log4 T L„NA(D! 1 lxV ) 



log R H 



x V 



1 We use "it" to denote an estimate of the random variable x. 
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Under relatively mild hypotheses on the distribution of N A , the two sets {dp* } and {d D x 2 } should con- 
verge to the same do when the mesh of the lattice goes to zero, that is, when R HxV — > oo. The rates of 
convergence, on the other hand, may of course differ. 
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Figure 3: Probability density functions (pdf ) for d£ */ for H x V = 33 x 32, 65 x 64, 1 29 x 1 28 and 257x 256. The 
dashed curve represents the pdf of the normal distribution with same mean and variance as the distribution 

for d 257x256 



*D,1 



The probability distribution function (pdf) for d^^ for critical percolation is shown in Figure|3]for four 
sizes of the lattice. The narrower distributions correspond to larger lattices (smaller mesh). The dashed 
curve is the normal distribution that has the same mean and variance as the distribution for d^ 57 ^ 256 . Up 
to an affine transformation, the distribution of logN A is the same as that of d{^ xV . If the lattice mo del a t 
(3 = 1 indeed has SLF_6 as its scaling limit, the Hausdorff dimension should be | by Beffara's theorem 1 2.2 i. 



The figure does not exclude this, but the convergence toward the Dirac mass 6 ( do — 2 ) is rather slow. Even 
for H x V = 257 x 256, log N A / log R is likely to take any value between 1 .6 and 1 .8. The standard deviation 

is therefore large for d^ xV and consequently for log N A . Note finally that, even though the pdf for d^ x v is 
skewed to the right, the window [d D * — 2s, d D * + 2s], where s is the standard deviation of the sample, 
is a fair approximation of the 95%-confidence interval. To determine the confidence interval for x 2 , we 
first compute it for N A and then take the log of both extremities of the interval. Taking the log should make 
the interval asymmetrical with respect to the mean. However, the intervals are small and this asymmetry 
can be ignored when estimating the confidence interval. 

Table [3] contains information on the confidence intervals for critical percolation CM{2,3) and the log- 
arithmic tricritical Ising model CA4{4,5). As discussed in Appendix |Bj critical percolation is the fastest 
to simulate and we have large samples. The logarithmic tricritical Ising model, on the other hand, is the 
slowest that we have studied and samples are considerably smaller. One can see that the confidence in- 
tervals for cIq^ and djj}^ are essentially identical. Note that the precision of measurements are very 
good with five or six significant digits. For our purposes, however, they are only acceptable. Indeed, sup- 
pose again that Beffara's result applies to the thermodynamic limit of these lattice models. Then, around 
H x V = 257 x 256, dp* remains more than 0.03 away from do = ? for £vW(2,3) (resp. 0.07 away from 
do = for CM (4, 5)), but between H x V = 257 x 256 and 513 x 512, this gap has been reduced only by 
0.002 (resp. 0.00032). The large variance and the slow progress of the mean towards the asymptotic value of 
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do make this measurement difficult. As the authors of [13] point out: "The asymptotic power law dependence 
of the various masses on R is approached relatively slowly." 

This makes the extrapolation H, V — > oo not only necessary but also delicate. How we perform this 
extrapolation is the topic of Paragraph|3.2.2| 
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___ 

d HxV 
U D,2 


|S| 


___ 

d HxV 
U D,1 


^ 

d HxV 
U D,2 


3 


2 


50000000 


1.64071|08 


1.69283|07 


100000000 


1 .52870106 


1 .58952|06 


5 


4 


1000000 


1.68060|36 


1.71113133 


100000000 


1 .58394|04 


1.61629|04 


9 


8 


1000000 


1 .69777\27 


1.71982|25 


50000000 


1 .60240103 


1.62413|04 


17 


16 


1000000 


1 .70754|21 


1.72500|20 


10000000 


1 .60952106 


1.62595|06 


33 


32 


1000000 


1.71396|18 


1.72856|17 


20000000 


1.61291|04 


1.62614|04 


65 


64 


1000000 


1.71915|15 


1.73164|15 


2201000 


1.61489|10 


1.62598|10 


129 


128 


1000000 


1.72312|13 


1.73399|12 


206000 


1.61609|27 


1.62564|27 


257 


256 


1000000 


1.72624|12 


1.73588|11 


229700 


1.61732|22 


1.62574|23 


513 


512 


1000000 


1.72867|11 


1.73731|10 


64400 


1.61764|38 


1.62513|38 


1025 


1024 


1000000 


1 .73074|09 


1.73855|09 









Table 3: d^}*, and d£j x 2 with confidence intervals for two of the models. The notation 1 .64071 108 means 
that the 95%-confidence interval is 1 .64071 ± 0.00008. 



3.2.2 Extracting d D from {d^J,, } and {dg* 2 v } 

By definition of the fractal dimension, the number of quarter-circles N a visited by the defect is proprotional 
to R d ° . Therefore, for sufficiently large R, one may assume that 

log R log R 



where a is another constant to be determined from the data. The approximation < |3.7[ | is linear in (log R) 
and a simple linear fit should be sufficient. However, the curves in Figure [2] clearly show that the sizes of 
the lattices we work with are too small for | |3.7) to be a reasonable fit. 

A better approximation is proposed in [30J where a renormalisation group analysis leads to 

logX ~ d ° + logR + R^logR (3 ' 8) 
where a, b are constants to be fitted and is a known exponent depending on the model. It is | for the 



Ising model, for example. Using 1 3.8 1 adds a hypothesis on the scaling behaviour of the model, that of the 
renormalisation group. Even though we did explore the quality of the fit obtained with | |3.8) , we use the 
following linear fits to get do : 



k 



log Na _ v- 1 



l ogR L- ^logR) 1 
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where the parameter Qo is to be interpreted as do ■ The number k of terms in the fit is chosen somewhat 
heuristically by requiring that the p-valu^]for each Qt remains below 0.05 and that k be the same for all 

models. We find that k = 2 is best for our data. 

Let us summarize some typical p-values obtained for fits of dp XV when the fits are based on all lattices 
but the two smallest ones: 3x2 and 5 x 4. In fits with functions 1 and j^-r, the p-values of the coefficient 

of 1 is always smaller than 10~ 14 and that of the coefficient of j^-^ is smaller than 4 x 10~ 4 , except in one 
case in which it is s» 0.3. When three functions are used, 1 , j^-^ and ( log R )2 / the p-values of the coefficients 
remain small: less than 10 _1 1 for the coefficient of 1 and less than 0.05 for that of r-^s- For that of ,, \^ 2 , 

logR (logR) 2 ' 

the p-value is also very small < 0.02 for all fits except those for the Ising model. That means that, for the 
logarithmic Ising model, the quadratic term might not be necessary. Indeed, in Figure |2j the curves for the 
logarithmic Ising model are seen to be exactly those where the concavity of the fit flips from being upward 
to downward implying that the coefficient of the quadratic term must be very close to zero. Finally, fits 
including a fourth term, ^p^p , have p-values larger than 0.1 for many of their coefficients, not only for the 

added term. 

We reject the values dj 4 ,*^ obtained for the smallest lattices. To choose which ones to reject, we proceed 
as follows. We reject the lattices whose inclusion (or rejection) cause a major change in the quality of the 
fit. Keeping the lattice H x V = 9 x 8, but rejecting the smaller ones, seems to be optimal. An indirect 
confirmation of this choice comes from the hypothesis that Beffara's theorem is valid for our lattice in the 
case of critical percolation. We also have to choose R. In the definition of the Minkowski dimension, R is 
the linear size of the box superimposed on the defect. An initial choice could therefore be the square root of 
the total number of quarter-circles: \/2H x V. However, the defect has access also to the half -circles on the 
boundary. A better choice is therefore the square root of the number of quarter-circles that can be occupied 
by the defect. In the present geometry, this is 



R = V2H x V + 4V (3.10) 

and this is the definition we use. This choice should not have a significant impact if large lattices are 
considered. For the sizes considere d here , though, there is an effect, as discussed below. The values reported 
in Tableware those obtained with (3.10 1 and k = 2 for both sets {d} 4 ,*^} and {d^^ v } based on lattices with 



H x V > 9 x 8 only. Using \3.8\ for the logarithmic Ising model gives do,i = 1 .6693 and do, 2 = 1 -6652; the 
quality of these results is comparable to the values in the table, but no better. As we shall see shortly, the 

limitation is clearly the errors on the d£5* V for large H x V and the choice of definition for the parameter R. 

We have explored also the possibility of predicting, from our data, both the fractal dimension do and 
the exponent of the correction term proposed by < |3.8) . If the correction term is small, this amounts to a 
nonlinear fit for the parameters en in 

N A =R a ° ai + (3.11) 

But we have only 1 lattices (1 1 for critical percolation) for 4 parameters. The Levenberg-Marquardt method 
(see [31 1, for example) does provide an excellent fit (in the sense that the curve is very close to the data), but 
the exponent do and 9 are totally off. Asikainen et al say that they use bins [Rt, Ri+i] with R^+i = \/2Ri- 
This gives them more than twenty data points. Still, like us, they claim that the prediction of the correction 
exponents is impossible with their range of data. We abandoned this idea. 

To evaluate the statistical error on the reported values doj and do ,2/ we perform the following simu- 
lation. We add noise to each datum of either set {dj 4 ,^ }, i = 1 , 2, distributed normally with zero mean and 
variance equal to the standard deviation measured from the previous experiment (see Figure [3}. With this 
new set of "noisy" data, we work out the fit as described above. Repeating this experiment a large number 



2 In a linear regression as the present ones, the p-value associated with one of the coefficients in the fit tests the hypothesis that this 
coefficient is zero. Therefore, a p-value of 0.05 for the coefficient c means that one might erroneously reject the hypothesis that c is 
zero in 5% of the cases. 
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of times (1000) gives a good idea of the deviation on dn.i or do, 2- For the best data, i.e. those for critical 
percolation, the noise coming from the width of the probability distribution of Na amounts to an error 

(= 95%-confidence interval) of 0.0006 on either dn,i or do, 2- This is the error quoted in Table |] We stress 
that this is the statistical error. But this method of evaluating the statistical error also gives an impression 

of the sensitivity of the measurements do.i and do, 2 on a single datum. We have repeated this experiment 
by changing one datum, that on the largest lattice H x V = 1 025 x 1 024, to the value at either extremity 

of its confidence interval. This changes do,i by a little more than two units of its fifth significant digit. In 
other words, the normal statistical error on a single datum leads to an error about a fourth of that on the 
whole set. This justifies our claim that the errors on the data for d{^ v significantly limit the precision on 
the dn.i's. Of course, it would be nice to probe lattices whose linear size is a few order of magnitude larger 
than the ones allowed by our slow MC algorithms. 



lattices used 


2 to 16 


4 to 32 


8 to 64 


16 to 128 


32 to 256 


64 to 512 


128 to 1024 


do ,2 with R 
do, 2 with R' 


1 .73484 
1.74210 


1 .74220 
1 .74430 


1.75013 
1 .74961 


1.75433 
1 .75348 


1.75113 
1.75115 


1 .74988 
1 .75003 


1 .74999 
1 .75004 



HxV 



Table 4: Measurements do, 2 for critical percolation CM{2, 3) obtained from a subset of the d D 

Janke and Schakel [14] report an interesting feature of their data for the Ising model. Their measurement 
of the dimension of the FK hull, the one to be compared with our dimension of the defect, shows the 
following behaviour. When they limit their analysis to their smaller lattices (L = 8 to 48), the dimension 

obtained is 1.665, very close to the predicted d^ 1 " 8 = I. But when they use their largest one (L = 64 to 
512), the result seems to slip down a little (to 1 .641). This is obviously a small effect. It should be noted 
that they report results, for this mass, only with the average cluster size x- In Table |3J we report the values 

do for critical percolation obtained from the set {dQ^ v } by extrapolation from 4 consecutive lattices. For 
example, the first fit is that with H x V = 3 x 2 to 17 x 16. (We admit, o f cou rse, that using 4 data for a 



3-parameter fit is not ideal.) The results of the upper line use the definition (3.10 of R as \/2H x V + 4V, the 
definition used throughout this analysis. These first move up and even overshoot the theoretical value, but 
then come back to extremely good ones. The problem reported by Janke and Schakel is not seen here. As 



underlined in Section 3.1 our settings are different. Our defect always crosses diagonally the lattice and one 
could guess that its gyration radius is fairly uniform. It remains to be seen how this could be an advantage. 
Note that we perform their experiment with our data for critical percolation as these go to V = 1 024. The 
same experiment done for our data for the logarithmic Ising model show the same behaviour (including the 
overshot and the final improvement of estimates when using larger lattices). Table|4]is completed by results 
done with respect to another definition of R, namely R ' = \/2H x V. Although this change hardly alters the 
fits on large lattices, it has an important effect on fits with small ones. This is to be expected. It nevertheless 
demonstrates how sensitive the extrapolation is to the definition of R. For the fits reported in Table |2j the 
switch from R to R' leads to a change in do within the confidence interval for all models but the logarithmic 
tricritical Ising model £A4(4, 5) where it is responsible for an increase twice as large as this interval. This 
indicates that, for the lattice sizes studied, there is no use trying to extract better extrapolations from our 
data, unless one can choose an R on a more theoretical basis. 
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3.2.3 Importance of boundary half-circles 



In ||32| , the role of boundary conditions as described by defects is studied for critical dense polymers 
£.M(1,2). The half -circles that are added along the boundary of the domain obviously change the prob- 
ability of the configurations. For example, if (3 > 1, configurations with many loops are favoured and 
those with several (shorter) paths joining two points of the boundary will be common. However, for crit- 
ical percolation that has (3 = 1, one may think that these half -circles do not play any role. This is not so. 
Configurations with and without these half -circles will have noticeably different do ■ 

Consider a long cylinder whose circular section contains an odd number of boxes. There will be (at 
least) one defect going from one end of the cylinder to the other. The boundary of the cylinder is now the 
union of two circles. The paths reaching these boundary circles can all be closed by boundary half-circles, 
all but one. Or they can be left open. In configurations without half -circles, the loops starting and ending at 
the boundary will be as many open loops competing for space near the boundary. The defect that crosses 
the cylinder will in this case be less dense close to the boundary. If, however, half-circles are added, they 
can be used by the defect to fill the space close to the boundary. Consequently, the dimension do of the 
defect measured on configurations without half -circles will be smaller than that measured on configurations 
with them. The latter is probably independent of the geometry. The former is likely to depend on it; for 
example, do depends on the length of the cylinder (without half-circles) and tends to do (measured with 
half -circles) when the ratio length / circumference tends to infinity. This interpretation is consistent with 
the growth process that is described by SLE. When an SLEg interface is grown on, say, the half-plane, it 
is allowed to bounce on the real axis (the boundary). The loop gas configurations with half -circles at the 
boundary are those allowing this reflection of the defect on the boundary. 

( a of intersections of defect 
with one line of boxes ) / (2H) 



0.50- 



0.45- 



0.40- 



0.35- 




Figure 4: Density of the defects on six cylinders as a function of the position along the cylinder axis. Bound- 
ary conditions are without (dashed curves) or with (full curves) half-circles. 

This effect of half-circles at the boundary can be seen in simulations even on long cylinders. In Figure |4j 
we plot the average of the ratio 

# of quarter-circles visited by the defect ^ 
total # of quarter-circles 

for each row of boxes along the cylinder. This ratio drops sharply near the boundary of the cylinder when 
half -circles are absent (dashed curves on the graph). When they are present (straight curves), the ratio is 
essentially constant on the full length of the cylinder. (The sizes given on the figure are labelled H x V where 
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H stands for the number of boxes along a section and V for the length of the cylinder.) The two ratios (with 
or without half-circles) differ significantly over a distance from the extremities that is about the length of 
the circumference. 

4 Other geometric objects: hull, cluster mass, external perimeter and 
red bonds 

4.1 Definition of masses for loop gas clusters 

The classical definitions of the various masses (cluster mass, hull, external perimeter and red bonds) need 
to be adapted to the context of a loop gas whose thermodynamical limit for rational k is argued to be a 
logarithmic minimal model. This extension is necessary as the random variables are not spins but boxes 
with only two states. It may provide a larger domain of exploration as the loop gas is parametrized by the 
continuous parameter k and not only an integer, like the Q -Potts models. 

Figure [5] (a) gives a configuration of the Ising model on a rectangular lattice. The reason for drawing it 
obliquely will become clear below. A configuration of FK clusters is obtained from this spin configuration 
as follows. Bonds between spins of distinct signs are erased and bonds between identical spins are removed 
(resp . kept) with probability ( 1 — p ) (resp . p ) where p = 1 — e ~~ a for a the inverse of the temperature and the 
Ising Hamiltonian of the form H = ^ Y.{ij)^<yi<yj — 1)- Figure|5j(b) shows a configuration obtained this way 
from the spin configuration in Figure [5] (a). The first step of a Swendsen-Wang upgrade consists precisely 
in constructing an FK configuration from a spin configuration. An FK configuration is partitioned uniquely 
into FK clusters, the maximal connected components of bonds and sites. We now recall the definitions of 
the various masses for these clusters. We shall concentrate on the largest cluster, the one shown in dark 
in Figure [5] (c). To distinguish between the classical definitions of these masses and the new ones to be 
introduced here, we add the label "spin" to the former and "loop" to the latter. 

The spin cluster mass is the number of sites connected by the bonds of the FK cluster. A site from the 
cluster belongs to the spin hull if a continuous curve can be drawn from this site to infinity (the boundary 
of the lattice) without crossing any of the cluster bonds. The spin hull is marked by pale (yellow) squares 
in Figure[5](d) and the other vertices of the cluster by dark (blue) squares. The spin hull mass is the number 
of sites in the spin hull, that is, the number of yellow squares. The hull partitions the whole lattice into 
three disjoint sets of sites: the hull itself, its interior and its exterior. Its interior consists of the sites of 
the lattice not belonging to the hull and from which there are no continuous paths to infinity avoiding all 
cluster bonds. The interior of the hull contains sites belonging to the FK cluster (drawn in dark (blue)) and 
others that do not. The exterior of a cluster is the complement of the union of the hull and its interior. The 
spin external perimeter (EP) is the subset of the exterior sites having at least one nearest neighbor in the hull. 
In Figure [5] (d), they are marked by dark circles (in green). The spin EP mass is the number of sites in the 
external perimeter. 

The last set to be described is that of spin red bonds (RB). Among the sites of the FK cluster, identify those 
that have the largest vertical coordinate and, among these, call cfnw the site with the smallest horizontal 
coordinate. Similarly, define as Gse the site with smallest vertical coordinate and, if there is more than one 
with this coordinate, the one among these with the largest horizontal coordinate. A bond in the FK cluster 
is a red bond if its removal breaks the FK cluster into two disconnected parts, each containing one of the 
two sites ctnw and crg^. In other words, if a voltage difference is applied between <7nw and crgE, the breaking 
of a single red bond will interrupt the electric current. The concept of a "red bond" is credited to Stanley. 
The definition appears in his paper 1 33 1 but not the name. If bonds are imagined as fuses, the first that will 
turn red (and burn) are the red bonds, hence the name. Note that the name "simply-connected bonds" is 
also used for the red bonds. This name is somewhat misleading as it does not have its usual topological 
meaning. It merely refers to the fact that these bonds keep the cluster connected. In the extension to loop 
red bonds, we shall apply the voltage difference between all topmost sites and all bottommost ones between 
the entry and exit points of the defects. This is a slight departure from the original definition that we have 
just introduced. The red bonds, obtained according to this new definition, are represented on Figure [5] (d) 
by small (red) slashes through them. 
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(a) (b) 




(c) (d) 



Figure 5: Four steps to construct various spin masses starting from a spin configuration. 



The correspondence between an FK configuration (obtained from a spin configuration) and a loop gas 
configuration is one-to-one, once the boundary condition of the latter has been chosen. We superimpose a 
grid whose mesh is 1 /a/2 times that of the spin lattice, as in Figure [6] (a). Spins of the first lattice occupy 
about half of the vertices of the new one. The state of each box in the loop gas configuration is determined 
from the FK configuration as follows. Two vertices of each box are occupied by spins. If there is a bond 
between these two sites in the FK configuration, the state of the box is such that the two quarter-circles do 
not intersect this bond. If there is no bound, then the two quarter-circles are drawn as if to prevent a bond 
to be drawn there. In such a loop configuration, the FK clusters not reaching the boundary are enclosed by 
loops as "tightly" as possible. 
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To construct the loop gas configuration corresponding to the FK configuration of Figure[5](b), we choose 
a boundary condition where exactly two defects enter at the top of the lattice and exit at the bottom. All 
remaining boundary boxes are joined by half-circles between neighbouring boxes. The result is shown in 
Figure [6] (b). Again, we define the loop masses for the largest object (given by the defects), even though our 
definitions hold also for the closed loops themselves. 

The geometric objects of the loop configuration, that we now propose to replace the classical (spin) 
masses by, are based only on the loop configuration and do not entail a choice of a compatible spin con- 
figuration (as one would in the second and last step of a Swendsen-Wang upgrade). They will not agree 
numerically with the spin masses of the original finite lattice but, hopefully, will have the same thermody- 
namical limit. 

The loop hull mass of any closed loop is the number of quarter-circles of the loop, including the two 
quarter-circles of each boundary half -circle that the loop visits. The loop hull of a defect is defined similarly 
as the number of quarter-circles from its entry to its exit. The loop hull created by the two defects is shown 
in pale (yellow) in Figure [6] (c). 

It is somewhat unfortunate that the word hull is used here as well as in the description of the SLE process, 
but with two different meanings. In SLE, the curve y, which is the counterpart of a defect, is drawn, say, 
on the upper half-plane EL There is no underlying lattice. The curve y is allowed to touch the real line (the 
boundary of EL) and itself. The hull drawn by y is the union of the points visited by y with the set of all 
points that have been cut off from infinity (in EL) by the curve touching the real line or itself. The SLE hull 
is therefore likely to include open domains of various sizes. In the loop gas, the defect never touches itself 
as it is the union of disjoint quarter-circles. The best approximation of a tangency point occurs in boxes 
whose two quarter-circles are visited by the same loop (or defect). Between the two visits of such a box, the 
loop has circled a region in the lattice whose interior could be naturally identified with an open set. The 
quarter-circles of this open set do not belong to the loop hull mass. This is contrasted by their counterparts 
in the SLE setting which do belong to the SLE hull. 

In the simulations, we will add the contributions of both defects, except for critical dense polymers 
with P = 0. Counting quarter-circles of both defects amounts to visualizing these two defects as left and 
right boundaries of a percolating cluster. For polymers, configurations with closed loops are forbidden. 
This means that the defects must visit every single quarter-circle of the lattice. Counting quarter-circles of 
both defects would then give the total number of quarter-circles for all (allowable) configurations and the 
fractal dimension of the loop hull mass would be exactly the predicted value, 2, and there is no need for 
simulations. Instead, we choose to restrict the counting to one of the two defects for this model (|3 = 0). 
This raises the question of whether the fractal dimension of defects depends on the number of macroscopic 
objects in the configuration, for example the number of defects crossing from top to bottom. This will be 
discussed in Section l4.3.3l 

The definition of the loop external perimeter (EP) is inspired by Grossman and Aharony's biased walker J34J. 
When a loop intersects an edge, it does so perpendicularly and its direction at this intersection is therefore 
one of four choices: north, east, south and west. We start by constructing the loop EP of the left defect. At 
the top entry point of this defect, the direction is southward. Consider the box the loop is entering. If only 
one of the two quarter-circles of this box is visited by the loop (the defect), the external perimeter follows 
this quarter-circle. If both quarter-circles are eventually visited by the loop, the external perimeter will draw 
a quarter-circle leading to the edge on its left whether or not this quarter-circle belongs to the state of the 
box. The loop has now progressed to a new face visited by the loop. Its next step is chosen in the same left- 
biased way. The loop EP mass is defined as the number of quarter-circles visited by the left-biased walker 
following the left defect from its entry to its exit. As we are going to explore two other ways of defining the 
EP, we denote the one just described by EP1. An example of a hull is drawn in Figure [7] (a) and its external 
perimeter, according to definition EP1, appears in Figure [7] (b). Note that the hull of the left defect is part 
of or to the left of its external perimeter. As can be seen there, the external perimeter cuts out many of the 
inner meanderings of the hull. To use Janke and Schakel's words [15 1, it is a smoother version of the hull. 
Even though we have not measured the external perimeter of closed loops, we complete our definition EP1 
for them. Find first the rightmost box edge crossed by the loop. If there is more than one such box, any 
of them will do. This edge is always horizontal. Starting from this edge with a southward direction, draw 
the EP by following a ngW-biased walker. Again, the hull will be part of or in the interior of the external 
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Figure 6: Four steps to construct loop masses. 



perimeter. The mass of the loop EP is again the number of its quarter-circles. It is smaller (or equal) to the 
mass of the loop hull. 

The next way of defining the external perimeter is closer to the interpretation of the two defects as 
enclosing a macroscopic FK cluster. In this interpretation both defects need to be considered and it is more 
natural to draw the EP for the left defect starting at the top and using a right-biased rule and then draw the 
EP for the right defect starting at its bottom and using again a right-biased rule. Figure|7](c) gives the result 
for the hull in Figure [7] (a). Note that it is also the external perimeter given as an example in Figure [6] (d) 
in mid-gray (green) superimposed on the hull in pale gray (yellow). As we shall see, this new definition 
(EP2) works equally well for all models but one. The problematic model is again critical dense polymers, 
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(c) (d) 



Figure 7: Three candidates for the loop external perimeter. Figure (a) shows the hull. 



the "extreme model" with |3 = 0. As said before, the defects of this model fill the lattice. Using a right- 
biased rule to obtain the EP of the left defect would draw an EP wiggling along the boundary, using only 
the quarter-circles touching the boundary. Drawing the EP of the right defect starting from the bottom and 
using a right-biased rule would have the same effect. This would give the same number of quarter-circles 
for all configurations, and this number grows linearly with the (linear) size of the lattice. The resulting 
fractal dimension for polymers would be 1 . This is not the predicted value and this definition does not 
probe the fractal nature of the geometric object. 

Our final proposal for defining the loop EP is denoted by EP3 and consists in adding the "inner bubbles". 
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The spin EP naturally adds sites that can be "deeply enclosed" within the cluster under consideration. 
Figure [5] (d) shows some of these bubbles. Indeed, besides the dark (green) disks to the left and to the 
right of the cluster, there is a handful of these disks trapped inside the cluster. We explore with EP3 the 
possibility of adding their counterparts in the loop gas case. Their addition requires several sweeps of the 
left- (or right-)biased algorithm. For example, after having drawn the external perimeter with, say, the 
right-biased definition for the left defect, we look for boxes visited by the hull but not by the EP just drawn. 
From those, we retain only boxes where one of the quarter-circles is not visited by the hull and where this 
quarter-circle is connected to the left lattice boundary by a path not intersecting the hull. Starting from 
this quarter-circle, we draw the right-biased EP of this "inner bubble". Figure|7](d) shows the resulting EP, 



with three "inner bubbles". As we shall see in Section 4.3.2 this last definition EP3 does not converge to the 



predicted value d^p but, instead, seems to converge to cIh- 

A loop gas configuration is a choice of state for each box of the lattice and these states draw loops, 
not clusters. To define clusters from these loops and subsequently a loop cluster mass, we first define an 
open bond. The original spin lattice has spins on only (about) half of the vertices of the loop lattice. We 
say that a box of a loop configuration has (or draws) an open bond if its quarter-circles do not intersect 
the segment joining the two vertices of the original lattice. Maximally connected sets of bonds are clusters 
and the loop cluster mass of such a component is the number of open bonds. This is similar to clusters in 
bond percolation from which the name "open bond" is borrowed. When the boundary condition allows for 
defects, we define the cluster between two (vertical) defects as the connected set of bonds between the two 
defects reaching both the top and bottom boundaries. This set is unique unless there is a closed loop located 
between the two defects and connecting both the top and bottom of the lattice. In the unlikely occurence 
of such a loop, we would use the leftmost connected set of bonds. The cluster mass is the number of its 
bonds. The geometric object used to define the loop cluster mass coincides with that of the spin cluster 
mass. However, the loop cluster mass counts bonds instead of vertices. 

Cutting a bond in the cluster inside a hull (as just constructed) amounts to flipping the state of the 
corresponding box. If one wants to break the cluster between two defects into two connected parts in such 
a way that the electric current from the top to the bottom is cut (as does the cutting of red bonds), this 
flipping of the box will have to change the two vertical defects into two defects, the first connecting the two 
entries at the top boundary, the second the two exits at the bottom boundary. This occurs if and only if one 
of the quarter-circles of the box being flipped belongs to the left defect and the other to the right one. We 
therefore define the loop red bond mass (or RB mass) as the number of boxes visited by both defects. As will 
be seen from Figures |5](d) and Figure [6] (d), the two definitions used to draw the spin and loop red bonds 
almost coincide, the difference being of one extra loop red bond in this particular example. 

4.2 Measurements 

The first thing to choose is the appropriate definition of R. This linear size is related to the size e of boxes 
in the Minkowski definition \3.1\ and appears explicitly in our approximation 1 3.4 1. It counts how many 



geometric building blocks (quarter-circles, bonds, sites) the mass under consideration may occupy. The 
natural thing to do is therefore to use two distinct definitions of R, one for the hull and the external perimeter 

Ri = V2H x V + 2H + 2V - 4 (4.1) 

counting the number of quarter-circles, and one for the cluster and red bonds 

R 2 = VH x V (4.2) 

counting bonds in the lattice. 

Figure [6] describes the setting of our simulations. Two defects enter from the top and exit from the 
bottom. Square lattices H x V are used with H a power of 2. For all models, the smallest lattice has H = 8. 
The largest lattice for critical dense polymers has H = 1 28, for critical percolation H = 1 024, while for 
all other models, the largest lattice has H = 256. The left defect enters in and exits from column number 
H/4+ 1, while the right defect enters in and exits from column number 3H/4. This choice leaves a quarter of 
the boxes to the left of the left entry point, and a quarter to the right of the right entry point. The algorithm 
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CM{\,2) 
CM{3,5) 
CM(2,3) 
£M(3,4) 
£M(4,5) 
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2.000 
1.833 
1.750 
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10 



1 1 



1.250 
= 1 .300 

1.333 
= 1 .375 

1.400 



C1EP 

1.250 
1.280 
1.326 
1.373 
1.380 



d* eo 



2.062 



g = 1 .925 1 .946 



£1-1.896 1.912 
-!# = 1 .875 1 .892 



f|§ - 1 .869 1 .884 



f = 1.250 



14 
15 



0.933 



| = 0.750 



13 
24 



0.542 



— = 425 

40 



dRB 

1.258 
0.956 
0.768 
0.557 
0.432 



Table 5: Fractal dimensions of dn (hull), dEP (external perimeter), dc (cluster) and dsc (simply-connected 



bonds). For critical percolation CM [2, 3), better estimates are available in Section 4.3 



used ensures that the defect entering at the top left will exit at the top right or at the bottom left. The entry 
points account for 4 edges of the lattice boundary; all other edges are matched pairwise by half -circles with 
a nearest neighbour edge. We limit our samples to configurations with defects crossing from top to bottom, 
like the configuration shown in Figure © The size of these samples are of at least 2 x 1 5 for each lattice, 
usually much larger for smaller lattices. For all models but critical dense polymers, we use the definition 
EP2 of the external perimeter. For critical dense polymers, we use EP1. 

The extrapolation is done as follows. For the hull mass, we do a fit with only 1 and 1 / log R as functions 
and using all lattices but the smallest 8x8. Figure [8] shows that the two functions 1 and 1 / log R are indeed 
sufficient. The fact that the data for critical dense polymers follow a slope different from the slopes of the 
other models can be easily explained. Recall that, for critical dense polymers, the two defects together fill 
the whole lattice and that, consequently, we choose to count only the quarter-circles of the left defect for this 
model. If we would do the same thing for one of the other models, the average number of quarter-circles 
would, by symmetry, be approximately half of what we actually measure. Therefore, by using formula 
( |3.6) , we would obtain 

.new lQg(N A )/2 log 2 

dH = logR =dH ~logR (43) 

as the new measurement of dn- It turns out that, for all other models, the slope of the fit is smaller in 
absolute value than log 2. This definition d™ w would therefore lead to fits with slope of the same sign as 
that obtained for critical dense polymers. For the three other masses, d^p, dc and dRB, one observes a clear 
departure from linearity. As an example, the five curves for the red bond masses are drawn in Figure [9] We 
therefore add a term 1 / (log R) 2 in the fit. As for du, we reject the smallest lattice 8x8 and, for dc, also the 
second smallest 1 6 x 16. The results of the measurements, using < |3.6) , appear in Table|5] 

The definition of loop hull mass that we have proposed is obtained from the number of quarter-circles 
visited by the two defects crossing the lattice from top to bottom. It is thus equal to the fractal dimension 
do measured in Section 3.1 if the object is the union of both defects crossing from top to bottom. If we 
are willing to accept that, in the thermodynamical limit, these two defects have minimal effect on one 
another, then dn should indeed be equal to do • This is confirmed by our simulations, and the hull is seen 
to have the best numerical estimate of the four masses considered here. This point will be discussed again 
in Paragraph |4. 3. 3| The agreement is less spectacular for the three other masses: dtp, dc and dRB. For these 
three, there seems to be a systematic error: the estimates for dEP are always smaller than the predicted one, 
and those for dc and dRB always bigger. It is for dc that the discrepancies are most obvious, especially 
because the five models considered have their predicted values all grouped in the narrow interval [1 .86, 2]. 
For example, ct^ 1 " 8 is closer to the d* eo for critical percolation than the one for the logarithmic Ising model. 

22 




1024 256 



64 



16 



H 



OLM(l,2) 
□ LM{3,4) 



512 128 
X LM(3,5) 
+ LM(4,5) 



32 



O LM{2,3) 



Figure 8: Measurements of cIh for the five models. 



In Paragraph |4.3.1 we argue that these discrepancies stem mainly from the small lattices our slow upgrade 
algorithms allow us to consider. 



4.3 Technical issues 

4.3.1 Estimates of critical percolation masses from smaller data windows 

As said before, we choose to reject the same lattices and use the same set of fitting functions for all models 
to obtain a given mass. But bigger lattices are accessible for critical percolation and better estimates for this 
model can therefore be obtained than those reported in Table [5] As we have done in Paragraph |3.2.2| we 
compute estimates d of the various masses from a window of only 4 lattices. For the four masses dn, d^p, 
dc and cIrb/ we choose to restrict the fit functions to 1 and 1/logR. The results are contained in Table|6] 
Only estimates using | |3.6) are reported. The agreement with theoretical values is better than reported in 
Table |5j significantly so for all but dn- These improvements support the hypothesis that agreement with 
predicted values could be improved also for other models by using larger lattices and, therefore, that the 
proposed definitions for loop masses have the same thermodynamical limit as their spin counterparts. 

FigurellOlshows the data for the four masses together with the fits obtained from the four largest lattices. 



4.3.2 The definitions of the loop external perimeter 

In Section 



4.1 



we proposed three definitions of the loop external perimeter that appeared to be natural ex- 
tensions of its spin counterpart. The first, EP1, counts the number of quarter-circles visited by a left-biased 
walker starting its journey at the top of the left defect and ending when this defect exits the lattice. The 
second, EP2, uses the quarter-circles visited by a right-biased walker traveling both defects, the left one 
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Figure 9: Measurements of cIrb for the five models. 



lattices used 


^theo 


8 to 64 


16 to 128 


32 to 256 


64 to 512 


128 to 1024 


d E p 
dc~ 
dRB 


I = 1.750 
| - 1 .333 
&~ 1-896 
| = 0.750 


1 .7502 
1.365 
1.815 

0.6880 


1 .7500 
1.342 
1.854 

0.7176 


1.7501 
1.333 
1.875 

0.7342 


1 .7502 
1.330 
1.886 

0.7428 


1 .7506 

1.331 

1.891 

0.7477 



Table 6: Measurements Z for the four masses of critical percolation CM(2, 3) obtained from a subset of the 
data. 

from top to bottom and the right one from bottom to top. Finally, EP3 adds to the quarter-circles of EP1 the 
"inner bubbles", cf . the discussion above. 

In our simulations, the definitions EP1 and EP2 appear to have the same thermodynamical limit and 
converge to the predicted value, whereas definition EP3 seems to behave in the limit as the hull itself. Even 
though critical dense polymers has the densest hull among the models considered, it cannot be used to 
explore the three definitions of EP as some of them fail for this model. Instead, we present here data for 
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Figure 10: Masses for critical percolation £.M(2, 3). Squares have been used exceptionnally for the cluster 
mass data to distinguish them from those of the hull. Fits are obtained using only the four largest lattices 
(H = 128,256,512 and 1024). Data for dn and dEp are compiled with a different definition of R than those 
of do and cIrb ; this accounts for the different horizontal alignment. 



£.M(3,5). After critical dense polymers, it is the model with the densest hull. The other models behave 
similarly to CM {3, 5). As the densest, the hull of £.M(3,5) is the one most likely to feel the boundaries 



and the two defects may impact each other. Figure 1 1 shows how the three definitions behave for lattices 



from H = 8 to 256. Fits with a quadratic term are obviously necessary, at least for EP1 and EP3; we use 
a quadratic term in the three fits, rejecting the smallest lattice in the three cases (and the H = 16 for EP2). 
The definition EP2 explores the side of the defects closest to the boundary and one might have thought that 
the hull would hit the boundary often; this would then bring its external perimeter EP2 closer to be linear, 
lowering its fractal dimension towards 1 . This is not the case. Actually, for all lattice sizes, EP2 remains 
larger than EP1, simply because it adds the contribution of both defects. If the density of EP2 is somewhat 
decreased by the boundary, our data do not reveal this. 

Obviously, the thermodynamical limit of the definition EP3 is not the predicted one with fractal dimen- 
sion || . For CA4{3, 5) and all other models, the extrapolation of the data we have seems to bring its limit 
very close to du- Because of the definition of EP3, dEP3 < cIh- As it nevertheless seems to overshoot cIh a 
little bit, a possible resolution is that the definition EP3 actually converges to the hull|^] If this guess is true, 
then it would mean that, as the lattice grows, the inner bubbles account for a fixed ratio of the hull. This is 
surprising because previous experiments did not see this. Of course, an explanation could be that the loop 
and spin external perimeters are physically different objects. But what about the definitions EP1 and EP2 
which do reproduce the dimensions of the spin external perimeter? A second possible explanation is the 
following. The external perimeter EP3 is that of one of the largest, if not the largest, geometric objects in 



3 The fact that it overshoots is of no real concern considering the extreme extrapolation that had to be done! Note that the data 
^EP3 V are m ' ne ran g e [1 -25, 1 -52] and the extrapolated value is 1 .88. This is not for the faint-hearted ... 
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the lattice. It is conceivable that inner bubbles play a significant role only in the largest objects. But in other 
experiments, fits are done on objects of all sizes within the lattice. (See for example the comment at the end 
of Section IV A. of [13].) Could the fits of these experiments have minimized the contribution of the largest 
objects? Yet a third explanation could be that, in adding the inner bubbles, we give them too much weight. 
Indeed, by comparing Figure [5] (d) and Figure [7] (d), we see that, where there is a single "inner spin" of the 
spin EP, we count 12 quarter-circles in the loop EP. Resolving this question requires more simulations. 

d E p 

2. 




Figure 11: The three definitions of the loop external perimeter for the model CM (3, 5). 



4.3.3 &h and d^p for configurations with one and two defects 

We have just seen that, for the left defect, EPs drawn by a right-bias walker (EP2) or by a left-biased walker 
(EPl) seem to have the same thermodynamical limit. This rules out a significant effect of the boundary on 
that geometric object. What about the influence of the right defect on the left one? One could argue that the 
right defect acts, for the left defect, as boundary on its right side and therefore should have no more effect 
on it than the boundary. This seems to be the case. 

Two sets of measurements are taken: one set on configurations with two defects (as we have done 
earlier) and a second one with only the left defect. The right defect is actually left in order not to change the 
size of the lattice, but is forced to remain along the top, right and bottom boundaries. The smallest sample 
is 10 5 . On these two sets of configurations, we measure dn of the left defect and its dEp with a right- 
biased and a left-biased rule. Figure 12 draws measurements of du and d^p of the left defect for the model 
CM. (3, 5). Measurements with two defects are represented by dashed lines. In all cases, they lie under their 
counterparts measured with only one defect. This is natural because, in the latter configuration, the defect 
has more space to wander in and it takes advantage of this. The four curves for the external perimeter 
are well separated over H = 128. The two obtained from a left-biased rule are over the other two with a 
right-biased one. Recall that the right-biased rule forces the external perimeter of the left defect to collide 
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often with the left boundary of the lattice. Consequently, the external perimeters thus obtained tend to have 
fewer quarter-circles than the two drawn with the left-biased rule. So the order in which the curves appear 
in the figure was predictable. But what is new is that the two ways of measuring cIh appear to have the 
same thermodynamical limit, and similarly for the four ways of measuring d^p. This confirms our intuitive 
statement. 
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Figure 12: Measurements for the model CM [3, 5) of cIh and dnp of the left defect with and without the 
right defect. 



Acknowledgements 

This work is supported by the Australian Research Council and the Canadian Natural Sciences and Engi- 
neering Research Council. 



A Loop gas and logarithmic minimal models 

We recall here the definitions of the loop gas and logarithmic minimal models CM{v,v') on me lattice. 
Each elementary face or box of a finite rectangular lattice, with H horizontal boxes and V vertical ones, is 
restricted to be in one of the two configurations in Figure 13 with the Boltzmann weights 



sin(A — u) sinu . . „, 

CJl = ^ — , CT = . (A.l) 

sin A sin A 



The crossing parameter A 6 (0, 5] labels the given loop model with loop fugacity 



P = 2 cos A. 
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(A.2) 



The parameter u is the spectral parameter. For real values, it measures the spatial anisotropy of the lattice. 
The simulations in the present paper are all carried out at the isotropic point u = A/2. In this case, the two 
box states in Figure 13 are equally probable, as cr = cri . 



(a) state 1 (b) state 

Figure 13: The two configurations of an elementary face with weights cri = sm ^^ u ^ and <Jq 



smu 
sin A ' 



There are many ways of fixing the boundary conditions. For the present discussion, let us join all loop 
ends at the four boundary segments to a nearest neighbour belonging to the same boundary segment. A 
4x4 configuration is shown in Figure 14 The Boltzmann weight of a configuration is determined by three 




Figure 14: A configuration on a 4 x 4 lattice. 



numbers: the number #i of boxes in state 1 , the number #o =HxV-#i of boxes in state and the number 
i of closed loops. The statistical weight is then 



P(#o,#i,«) =dt'cr* (3VZ 



(A.3) 



where the normalization factor (the partition function) Z assures that the sum of the P(#o,#i ,£) over all 
2 Hx v configurations is 1 . In the isotropic case, the lattice model depends on the single parameter, the loop 
fugacity (3, and the probability of a configuration is determined by the number of loops that it contains. 
Care has to be taken if |3 = since the partition function itself vanishes with the boundary conditions 
indicated in Figure 14 This does not concern us here since the boundary conditions used in the bulk part 
of this paper all involve at least one so-called defect. 

As a CFT, the generic loop gas is described by the Coulomb gas Il23l . When ^ is rational, however, the 
model has some remarkable properties. Assuming < A < n, there then exist coprime integers p and p ' 
such that p ' > p and 

V 



p' K 



This choice of crossing parameter defines the logarithmic minimal model £A4(p,p')- It is argued in [24J 
that the continuum scaling limit of the loop gas with these special values of A yields a logarithmic CFT |25 1, 
denoted by CM. (p, p ') and with central charge 



1 6A 2 6(p'-p) 2 
c = 1 = 1 . 

7t(7T— A) pp' 
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(A.5) 



Particular examples of these continuum limits correspond to critical dense polymers CM. (1,2), critical per- 
colation CA4 (2,3) and the logarithmic Ising model CA4 (3,4). It is noted that, even though this paper reports 
measurements only for rational X/n, the definition of a loop gas on a finite lattice holds for every A. 



B Upgrade algorithms 
B.l Critical percolation ((3 = 1) 

Critical percolation is described by CM.{2, 3) where A = j and c = 0. For the purpose of simulations, this 
is the simplest case as (3 = 2 cos A = 1 and all configurations are equiprobable at the isotropic point where 
u = j — ^ . At this point, we sample the configuration space according to ( |A.3[ ) by choosing randomly the 
state of each box of the H x V lattice, giving each state in Figure [13] equal probability. Since |3 = 1, no loop 
counting is required and obtaining good samples on large lattices can be done swiftly. 

B.2 Models with < (3 < 2 

For A G [0, j )/ the loop fugacity |3 ranges in (0,2]. For (3 < 1 , the probability P (A3) can be seen as penalizing 
configurations with a large number of loops. On the contrary, for (3 > 1, configurations with many loops 
are favoured. For [3 ^ 1, an upgrade algorithm is needed to sample the space with the correct probability. 

We recall that the Swendsen-Wang algorithm takes advantage of the FK representation of the Potts 
model. For spin models, like the Q-Potts model, the Swendsen-Wang algorithm provides an efficient way 
to upgrade configurations. From a given spin configuration where each site is in one of the Q available 
states, the Swendsen-Wang algorithm first constructs an FK configuration by removing bonds between 
neighbouring sites that are not in the same state and by keeping, with probability p = p ( Q ), bonds between 
neighbouring sites in the same state. The FK graphs are in one-to-one correspondence with the loop gas 
configurations. The second step of the algorithm simply chooses one of the Q possible states for each of 
the connected components of the FK graph. But this is meaningful only if Q is an integer. The relationship 
between the number of states Q of the Potts model and the fugacity is Q = (3 2 . Except for (3 = 1, n/2, \ 3> and 
2, the corresponding Q of the loop gas will not be an integer and the second step of the Swendsen-Wang 
algorithm cannot be extended in an obvious fashion. 




(a) no entering defect (b) one entering defect (c) two entering defects 



Figure 15: The effect of flipping one box on the number of loops. 



We therefore use the traditional Metropolis algorithm, flipping one box or a small number of them at 
a time. Suppose that we have allowed defects to enter and exit the configurations. Flipping one box will 
always change the number of loops by 0, +1 or — 1 . This is demonstrated in Figure [15] The first pair shows 
the situation where the two quarter-circles of the box being flipped belong to closed loops. The portion of 
the closed loops outside the box may contain a large number of quarter-circles. We represent them as small 
"ears". The configuration before the box flip is one of the two configurations represented in Figure 15 (a) 
and the configuration after is necessarily the other one. The change in the number of loops is therefore ±1 . 
Figure 15 (b) depicts the case where one of the two quarter-circles belong to a defect. The external parts of 
the defect are represented as segments. Again, the change is ±1 loop. If two distinct defects cross the box as 
in Figure 15 (c), its flipping does not change the total number of loops. Note that if the same defect crosses 
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the box twice, the situation is that of the right drawing of Figure 15 (b) and a loop will be created by the 
flipping process. 

The simplest MC algorithm is therefore to flip one box at a time. The box is chosen at random in the 



H x V lattice. If the number of loops has not changed (as in Figure 15 (c)), we accept the flip. If the number 
of loops has changed, the ratio 

P(configurationafterflip) _ ± 1 
P (configuration before the flip) 

is computed. If this ratio is larger than one, the change is accepted. If it is smaller than one, the change is 
accepted with probability given by the ratio. If this algorithm is used repeatedly, the probability of visiting 
any configuration starting from any other is nonzero. This is the usual Metropolis algorithm. 

In order to calculate the change in the number of loops, it suffices to construct the loops intersecting the 
box being flipped. All our simulations involve one or two defects and these are large. It might therefore 
be useful to flip several boxes at a time, as several could intersect defects. There is a compromise to be 
struck between changing many boxes, to speed up the upgrading, and maximizing the acceptance rate. 
Fortunately, the acceptance rate does not fall too quickly as the number of boxes being flipped is increased. 
This can be seen as follows. Let x denote one of the H x V boxes of the lattice and let a x be a given 
configuration on all boxes but x. The choice of state for the box x gives rise to two configurations for 
the whole lattice. If both defects cross box x, then the acceptance rate is 1 since the number of loops is 
unchanged. Exclude now the case of two defects crossing x. The difference between the numbers of loops 
between these configurations is then ± 1 . If the configuration with the largest number of loops is cr J and 
the other cr~, the ratio of their probabilities is 

and the sum of the conditional probabilities is 1. We thus have 

p ( ff xl«-x) = p^y and p «l^=pTr (B ' 3) 

Suppose now that |3 < 1 and then P(c x |<r x ) > P(crJ|cr x ). If we start from crj, the flip of x will be accepted 

with probability 1 and, if we start from cr~, it will be accepted with probability (3. The acceptance rate is 
therefore 

r = P(a+|a x ) + (3P(a-|a x ) = -^ T , p<1. (B.4) 

Likewise for |3 > 1 , one gets 

r = 2/(1 + 0), p>1. (B.5) 

The acceptance rate is therefore 

r = 2min(P«|cT x ),P(cT-|(j x )). (B.6) 

What is the acceptance rate when several boxes are flipped at once? Since a loop could cross more than 
one of the boxes flipped, the analysis is more delicate. But if the number n of boxes flipped is much smaller 
than H x V, the hypothesis, that no loop crosses more than one of these n boxes, is reasonable. Let cr~ be 
the configuration with the smallest number of loops among the 2 n configurations under consideration. Let 
p _ be its conditional probability, that is, its probability given that the other ( H x V) — n boxes are held fixed. 
If (3 < 1 , then cr~ is the most probable and a calculation similar to the previous one leads to the acceptance 
rate of 

z (:>"-)- 

V ' 7 0<i<n/2 V 7 J 

The result is the same if (3 > 1 . In either case, the probability p_ is 



p_ = (1 + (3)- u . 
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As an example, the acceptance rate for the logarithmic Ising model ( (3 = y/2) of the flip of a single box is 
82% and that of the flip of 20 boxes is 45%, always under the assumption that no loop crosses more than 
one of these 20 boxes. On large lattices (H x V = 128 x 128 or larger), the measured acceptance rate f is 
the one calculated, within measurement error. The assumption on loops crossing the boxes being flipped 
is thus confirmed as being reasonable. By flipping several boxes at once, we save some overhead in the 
upgrade process and time when more than one of the n boxes intercept the defects. 



In Section 4.2 we study configurations with two defects. There are two entry points at the top of the 
lattice and two at the bottom. Our study conditions the defects to enter from the top and exit at the bottom. 
The algorithm is not forcing the production of configurations with defects running vertically only. That 
is, the algorithm gives also configurations where the defect entering at the top left entry point exits at 
the top right one. We reject these and keep for measurements only those where the defect entering at 
the top left entry point exits at the bottom left exit point. Since upgrades are (time) expensive, we have 
explored the following way of gaining time. After the completion of the number of upgrades necessary for 
independence, we check whether the defects are vertical. If they are not, we perform upgrades (checking 
between each) until the new configuration has vertical defects. This unfortunately biases the MC process. 
Indeed, choosing the first configuration with vertical defects after a run of configurations without them is 
likely to shift the distribution, probably onto configurations that have a constriction point (where the two 
defects have switched from a left-right to top-bottom pattern). To avoid this bias, a configuration without 
vertical defects is simply rejected and a full number of upgrades is redone before checking for vertical 
defects and possible measurement. 

B.3 Critical dense polymers (|3 = 0) 

This case is particular because, for |3 = 0, loops of any size are forbidden. Defects are thus necessarily 
present and must fill the lattice. We have developed a (very slow) algorithm to generate randomly dense 
polymer configurations with two defects. The difficulty here is not to sample the configuration space with 
the correct probability. As for critical percolation, all (admissible) configurations are equiprobable. The 
problem is to make sure that they do not contain loops. Suppose that we are interested in configurations 
where two defects enter at the top edge of a rectangular lattice and exit at the bottom edge, and that a con- 
figuration without loops is already known. The entry and exit points are fixed throughout the simulations. 
The algorithm first lists all boxes that are crossed by both defects. This list is never empty. One of these 
boxes is then chosen randomly with uniform probability and flipped. The type of flip is the one depicted 



in Figure 15 (c) and the resulting configuration is also without loops. However, one of the defects enters 
from and exits at the top while the other goes bottom to bottom. By repeating the algorithm we get a new 
configuration with no loops and two defects crossing the lattice from top to bottom. 

How does one get a starting configuration without closed loops? We fill the lattice randomly. The result 
is a configuration with (in general) many loops and two defects. We then repeat the following step. We list 
boxes that are crossed only once by either defect. If there is a closed loop, this list is non empty. We then 



flip one of these boxes, thereby changing the left configuration of Figure 15 (b) to its right companion, and 
start over. At the end we can check that the defects cross from top to bottom and, if necessary, correct by 
(half) the above algorithm. 



C Initial thermalization and independence of measurements 

Let us call an MC step the trial of flipping 20 boxes at once. We have established experimentally the number 
of MC steps necessary to assure that measurements of fractal dimension of defects or of other masses are 
statistically independent. As is well-known, this number of steps strongly depends on the model. For 
example, the 4-Potts model requires a significantly larger number of steps than the Ising model. For large 
lattices, we use the following number of MC steps between measurements 
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where d(|3) is a factor dependent on the model: d = 1 for CM{3,5), d = 2 for the logarithmic Ising model 
CM (3, 4) and d = 4 for the logarithmic tricritical Ising model CM (4,5). Hence, for CM (3,5), we try to flip 
one eighth of the boxes (by trials of 20 boxes) between measurements. This number goes up to one half for 
CM(4, 5). A finer analysis is under way and will be published elsewhere. It confirms that these numbers 
are fine, probably on the zealous side. 

The initial random configuration is thermalized by 1 000 times the number of MC steps just listed. The 
initial thermalization is therefore four times longer for CM(4,5) than for .M(3,5). Again, we think that 
this is satisfactory, if not "too prudent", but this is harder to assess. We have experimented with the (spin) 
Ising model and compared it with the |3 = v2 loop gas as follows. First, we use an ordinary Ising spin 
lattice wrapped on a cylinder at the critical temperature. There are four times as many spins along the 
cylinder than along the circular sections. We measure the two-point spin correlation function on the section 
in the middle of the cylinder using many samples of 1 3 configurations (chosen to be independent). We 
then perform the same experiment with the corresponding loop gas. We know that a loop configuration 
can be interpreted as an FK configuration and therefore that a spin configuration can be constructed from 
it. We repeat several times the following experiment: from a random configuration, we apply N MC steps 
and then proceed to measure the spin correlation function with 1 3 configurations, each separated by the 
number of MC steps assuring statistical independence. By comparing the distributions of the correlation 
function for the spin model, on one hand, and the loop gas on the other hand, we find that it is sufficient to 
use the N stated earlier, that is N = 1 000, times the number of MC steps between measurements. 

For critical dense polymers ((3 =0), no initial thermalization is necessary. We flip H x V/8 boxes between 
measurements . 

As the boxes are chosen randomly (for critical dense polymers and the other models), the number of 
boxes actually flipped might be smaller since, for example, the same box can be flipped twice between two 
measurements. Moreover, if a box is chosen more than once among the 20 being flipped, it will be flipped 
only once and the total number of boxes flipped during this Metropolis step will be smaller than 20. This 
assures that the number of boxes flipped at each step can be of either parity. 

We use L'Ecuyer and Tezuka's pseudo-random number generator [35J in the simulations. 
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